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Singlet-triplet dispersion reveals additional frustration in the triangular dimer 

compound Ba3Mn208 
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We present single crystal inelastic neutron scattering measurements of the S — 1 dimerized 
quasi-two-dimensional antiferromagnet Ba3Mn208. The singlet-triplet dispersion reveals nearest- 
neighbor and next-nearest-neighbor ferromagnetic interactions between adjacent bilayers that com- 
pete against each other. Although the inter-bilayer exchange is comparable to the intra-bilayer 
exchange, this additional frustration reduces the effective coupling along the c-axis and leads to a 
quasi-two dimensional behavior. In addition, the obtained exchange values are able to reproduce 
the four critical fields in the phase diagram. 

PACS numbers: 75.10.Jm, 75.40.Gb, 75.30.Et 



The presence of exotic ground states [l[ , excitations [4| 
and quantum phase transitions [3| in a growing number 
of available materials is making low-dimensional gapped 
frustrated quantum magnets especially topical. In some 
cases, the ability to tune the system Hamiltonian via ex- 
ternal fields allows access to additional phases. For in- 
stance, the singlet-triplet gap in the spectrum can be 
closed via application of external magnetic fields [J, Q . 
The resulting quantum critical point (QCP) separates a 
low field disordered paramagnct from a high field ordered 
phase which to a good approximation can be considered 
a Bose- Einstein condensation of magnons [y, |7| ■ 

The hexagonal antiferromagnetic bilayer Ba3Mn2 08 
is established as a quasi two dimensional (2D) frus- 
trated antiferromagnet with a quantum critical phase 
diagramjsl, Isj. Rather than a single Bose-Einstein con- 
densate phase as in 5 = ^ singlet-triplet antiferromag- 
nets, Ba3Mn208 has two sequential magnetically ordered 
phases, i.e. four zero-temperature QCPs, as a function 
of applied magnetic field: ^qHci ~ 9 T, u,riHr.2 ~ 26 T, 
fioHcs ~ 32.3 T and ^o^c4 ~ 48 T[il[li,El. The first 
(second) ordered phase for Hd < H < Hc2 (ffc3 < H < 
Hc4) can be approximately described as singlet-triplet 
(triplet-quintet) condensation. The nature of these high 
field ordered phases fundamentally depends upon the 
zero field exchange paths and constants. 

Powder INS measurements of Ba3Mn2 08 yielded three 
independent exchange constants through comparison to a 
spherically averaged scattering intensity [10|. The dimer 
exchange Jq and interdimer exchanges are illustrated in 
Fig. [Ijb). Examining these excitations using single crys- 
tal INS, we find next nearest neighbor (NNN) interac- 
tions J4 play a significant and heretofore unexplored role 
in the propagation of triplet excitations in Ba3Mn208. 
Since the dispersion bandwidth only depends on the dif- 
ference J2 — J3, we use ^i^Hci 12,ll3illJl to determine the 
J3 exchange. We find that both inter-bilayer exchange 
constants, Ji and J4, are important in determining the 
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FIG. 1: (a) Constant Q scans in the (hhl) and (hkO) planes. 
Solid lines are Gaussian fits. The (|00) fit includes a fixed 
elastic Gaussian and a higher energy Gaussian to account for 
incoherent elastic scattering and a scattering angle dependent 
background respectively. Data are offset vertically for pre- 
sentation, (b) Crystal structure of Ba3Mn208 showing Mn^""" 
(red) and 12-fold coordinated Ba sites (green) [ill]. Oxygen 
sites (blue) involved in coordination of two dimers are shown 
to illustrate interdimer connectivity and tetrahedral coordina- 
tion of the Mn^^ . 10- fold coordinated Ba sites are not shown 
for clarity. Exchange connectivity within and between S — 1 
dimers is shown with different line types. J3 is only shown 
for a single pair of dimers. Exchange constants correspond to 
Mn-Mn distances of 3.985 A(Jo), 4.569 A(Ji), 5.711 A(J2), 
6.964 A(J3) and 7.313 A(J4). 



critical fields. In addition, we show that although fer- 
romagnetic, Ji and J4 are competing interactions due 
to the geometric frustration of the triangular bilayers 
and the inter-bilayer exchange paths. This explains why 
Ba3Mn2 08 is a quasi-2D magnet in spite of the fact that 



the dominant inter-bilayer coupling, Ji, is larger than 

\J2\- 

Single crystals of Ba3Mn208 were grown from a flux 
melt as described in Ref. |9|. Five crystals with a net 
mass 0.372 g were coaligned using the HBIA triple-axis 
spectrometer (TAS) at the High Flux Isotope Reactor 
with a resulting mosaic of « 42'. INS measurements 
were performed on the SPINS cold neutron TAS at the 
NIST Center for Neutron Research. The instrument 
was configured with 80' coUimation before the sample 
followed by a cooled Be filter and an 80' radial col- 
limator. Measurements were performed with fixed fi- 
nal energy Ef — 5 meV selected by a 10.5 cm wide 
horizontally focussing analyzer. This configuration re- 
sults in a measured elastic full width at half maximum 
(FWHM) energy resolution at the (0.3 0.3 3) wave vec- 
tor of STiLu = 0.294(7) meV. The majority of our data 
presented here consists of constant wave vector Q scans 
performed at T = 1.7 K in the {hhl) scattering plane. 

Representative constant Q scans are shown in 
Fig. [TJa) . The data show a single resolution limited mode 
propagating along the (/i/i|) direction between energy 
transfers hco of 1 and 3 meV. The bandwidth and gap 
of the spectrum is in good agreement with values deter- 
mined from powder INS measurements [10|. Measure- 
ments at several high-symmetry locations in the (hkO) 
scattering plane were also performed. These were made 
in the identical instrumental configuration with the fo- 
cusing analyzer broadened to w 23.1 cm to integrate over 
additional scattering intensity. The broadened lineshape 
for (^00) in Fig. [IJa) reflects the broadened Q resolution 
of the wider analyzer. As shown in Fig. [Ij^a), the mea- 
sured scattering intensity in the {hkO) plane is very weak 
compared to the (hhl) plane. 

To illustrate the dispersion and wave vector depen- 
dent scattering intensity, we combine a series of constant 
Q scans in the image shown in Fig. [2I^a)-(d). Figure [2{i) 
shows the path through wave vector space for these data. 
Figure [H indicates a single dispersive mode in both the 
{hhQ and {((l) directions. The dispersion along the 
(hhd^) wave vectors has a distinct W shape with clear 



minima in the vicinity of h 



i and h 



and a local 



maximum at the h = i^. The dispersion along the I wave 
vector has a periodicity of three reciprocal lattice units 
(rlu). 

126 individual constant Q scans were separately fit to a 
single inelastic Gaussian peak. The peak positions along 
several high-symmetry directions are plotted in Fig. [2ja)- 
(d) and along several (CCO wave vectors in Fig. [3l[a). 
The singlet-triplet excitation is fairly dispersive along the 
(OOZ) direction as shown in Fig. [2{d) and also along the 
/ direction near the zone boundary as shown in Fig. [3l^a) 
for (0.15 0.15 I). However, there is almost no dispersion 
along I for 0.3 < h < 0.7. In fact for the (CCO scans 
shown in Fig. EJa), the dispersion fiattens and continu- 
ously changes phase as one moves from C = to C = ^ 



The three original exchange constants determined from 
analysis of powder INS data do not account for the lack 
of dispersion along the I direction near the zone center. 
Rather we must include the NNN inter-dimer interaction 
in the Heisenberg Hamiltonian J4. As in the powder 
analysis, we use a dispersion based upon weakly coupled 
S = I dimers: 



nio^Q)^jAl + -A,J{Q)RiT) 



(1) 



where v — {1,0,-1} denotes the spin projection S^ of 
the triplet excitation, J'(Q) is the Fourier sum of in- 
teractions beyond dimer exchange, Aq = Jo + 2Z3/3, 
A±i — Jq — D/3 and R{T) is the thermal population 
difference between ground and excited states. For 5 = 1 
antiferromagnetic dimers. 



R{T) = 



1 - exp(- Jo/3) 

- 3 exp(— Jo/3) -I- 5 exp(— 3 Jo/3) ' 



(2) 



where f3 — j-^- Including the J4 interaction, the Fourier 
sum becomes 
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(3) 



The small anisotropy D = —0.032 meV [9|, llOlha s been 
determined from zero-field ESR experiments [l5|. For 
Ba3Mn2 08 in zero magnetic field, D will not change 
the dispersion and will only broaden the constant Q 
scans. As shown in Fig. [IJa), J4 corresponds to NNN 
interdimer interactions along the c-axis between bilay- 
ers. This exchange path is through the 12-fold coordi- 
nated Ba sites which are co-planer with the O sites as- 
sociated with the Mn tetrahedra. Each of the 12-fold 
coordinated Ba sites participates in three such J4 ex- 
change paths. A simultaneous fit of the peak positions of 
the constant Q scans in the {hhl) scattering plane yields 
Jo = 1.642(3), Ji = -0.118(2), (J2- J3) = 0.1136(7) and 
J4 = —0.037(2) meV. J2 and J3 represent interdimer ex- 
change between the same spin dimers. They cannot be 
independently determined from the zero-field dispersion 
because J'(Q) depends on their difference. The phase 
shift for triplets propagating between dimers via J2 or 
J3 (see Fig. [IJb)) is due to the opposite symmetries of 
the singlet (symmetric) and the triplet (antisymmetric) 
states under a permutation of the two sites that form 
the dimer. The corresponding dispersion curves shown in 
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FIG. 2: (a)-(d) T = 1.7 K scattering intensity of Ba3Mn20g versus hcj and Q. Figure is interpolated from constant Q scans 
in the (hhl) plane. White circles are peak positions from Gaussian fits described in the text. Solid line is the fitted dispersion. 
(e)-(h) Integrated scattering intensity from Gaussian fits to constant Q scans. Solid line is the fitted scattering intensity, (i) 
Path through the (hhl) plane shown in (a)-(h). Dashed lines correspond to scans shown in Fig. [S] Contour lines are the fitted 
dispersion with dark[light] lines representing smaller [larger] energy transfers shown at 0.2 meV intervals. Horizontal axes of 
(a)-(h) are plotted in reciprocal lattice units (rlu). 



Figs.ma)-(d) and [3]Ja) agree very well with the measure- 
ment. We also plot the fitted dispersion along (|^/) for 
J4 = (dashed line in Fig.JSlJa). It is clear that the NNN 
exchange is needed to describe the measured dispersion. 
The dispersion minimum hu = 1.081 meV occurs at wave 
vectors (i-(5, i-(5, 1) and (| + (5, j + 6, 1) with 5 « 0.0259. 
This value is similar to that determined from our powder 
measurements, and recent mean-field calculations [IJ]. 

Figures [2fe)-(h) and[3Ub) show the Q dependent inte- 
grated scattering intensity. Based upon the established 
SMA dimer modelflO], the magnetic scattering intensity 
for fcsT < Jois 

cos(Q ■ d) 



/,„(Q,w)cx|F(Q)|- 



1 



-5{hLu-nLu{Q)), (4) 



where \F{Q)\ is the magnetic form factor [16|,|l7| and d is 
the intradimer vector: d = 0.1858c. The dimer structure 
factor explains the lack of intensity found in the {hkO) 
plane, i.e. Q^hko) -L d. We fit the integrated intensity 
to the SMA model in Eq. (|4]) using the dispersion found 
earlier and only a single multiplicative fitting parame- 
ter. The fit agrees quite well with the measured results. 
The dimer structure factor accounts for almost all of the 
scattering intensity. 



The interdimer exchange corresponding to triplet prop- 
agation out of the hexagonal a6-plane is ferromagnetic 
for both Ji and J4, whereas interdimer coupling in the 
hexagonal a6-plane is entirely antiferromagnetic and con- 
sequently geometrically frustrated. We can gain insight 
into the consequences of the two interlayer couplings by 
considering the magnetic structure for Hd < H < Hc2- 
The in-plane geometric frustration alone leads to a 120° 
magnetic structure for a single isolated bilayer. Such a 
structure blocks the inter-layer couplings Ji and J4 be- 
cause a spin Si on a given bilayer interacts with the sum 
of three spins that are on the three different sublattices 
(A, B and C) of the adjacent bilayer (this sum is zero 
for a 120° structure). Consequently, the ferromagnetic 
interaction Ji induces a spiral phase {a ^ 120°) [9|, [l^] 
to have a net moment on each triangle. This moment 
is ferromagneneticaly coupled to the spin Si that is at 
the center of the corresponding triangle (and in the next 
bilayer). However, such spiral distortion is not favorable 
for the ferromagnetic J4 term: the three spins that are 
next-nearest-neighbors of Si in the next bilayer form a 
bigger triangle and the net moment of this triangle is 
antiferromagnetically aligned with S,. Therefore, the Ji 




to enhance the quasi-2D behavior of the bilayer. The 
refined exchange constants also provide a zero-field con- 
firmation of the high magnetic field QCPs. In addition to 
marked progress in understanding the quantum critical 
phase diagram of Ba3Mn2 08 and its excitations and 
ordered states, the determined dispersion should also 
be applicable to other isostructural gapped magnetic 
compounds in the A3Cr20s series where A =Ca^"'", Sr^+ 
or Ba2+ such as BaaCrzOg [H, [ll 0, [2l| and SrgCrzOs 
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FIG. 3: (a) Dispersion and (b) integrated scattering intensity 
of Ba3Mn2 08 determined from Gaussian fits to constant Q 
scans. Solid lines in (a) are fits to the dispersion as discussed 
in the text. Dashed line in (a) is the fitted dispersion for 
i^^l) fixing J4 = resulting in exchange constants Jo — 
1.640(8), Ji = -0.117(5) and J2 - J3 = 0.109(2). Solid lines 
in (b) are fits to the scattering intensity as described in the 
text. 



exchange interaction competes against J4 producing a 
flatter singlet-triplet dispersion along the c-axis, i.e., an 
effective quasi-2D like behavior. This competition also 
leads to smaller deviation of a from 120°. 

Based upon the determined exchange constants and 
an appropriate choice of J3 = 0.167 meV such that Hd 
agrees with its well determined experimental value, we 
are able to determine the remaining critical fields. Using 
a two-level approach I4I keeping only the singlet {S^ — 2 
quintet) and the S*^ = 1 triplet dimer states for Hd < 
H < Hc2 {Hc3 < H < Hci), we obtain values for the four 
critical fields: H^ = 10 T, Hc2 = 27 T; H^j. = 32.75 T; 
Hc4, = 47.75 T. The agreement with the experimental val- 
ues is quite good considering that neglecting the excited 
triplet S^ = 0,-1 states introduces errors of the order 
of 1 T. For instance, if we include the effect of quadratic 
fluctuations to the S*^ = 0,-1 triplet states, the lowest 
critical fleld becomesiJci = ^^iQmi7i)/gfJ'B = 9 T in ex- 
cellent agreement with the experimental value [Qmin is 
the wave- vector that minimizes w^(Q)]. 

Through single crystal INS measurements of 
Ba3Mn208 we find that NNN exchange interactions are 
significant in describing the system. The nearest and 
NNN inter-bilayer couplings are both ferromagnetic. It 
is clear that the antiferromagnetic interactions within 
the bilayer lead to a geometrically frustrated system. 
However for Ba3Mn208 we find a novel case where 
ferromagnetic inter-bilayer exchange actually serves 
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